Condensed states of a semiflexible copolymer in poor solvent: Figures of eight and 

discrete size torii 
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We examine the condensed states of a simple semiflexible copolymer in which there are two 
monomer types that are immiscible with each other and with the solvent. Although this is similar 
to the well known problem of collapse for a semiflexible homopolymer we find that it gives rise to 
a much richer variety of condensed states. We predict the existence of these states using simple 
analytic arguments and also observe them directly using Brownian dynamics simulations. 
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One of the most fascinating physical properties of DNA 
is its collapse from a diffuse coil to a compact toroidal 
structure under the addition of condensing agents to the 
solvent. This phenomenon is not only of biotechnologi- 
cal importance [l|, |3| but also serves to illustrate a gen- 
eral principle of competition between bending and sol- 
vent energy that is common to the collapse of all semi- 
flexible polymers including other biopolymers, synthetic 
polymers and also carbon nanotubes p] . Although much 
progress has already been made in understanding the ki- 
netics of folding IJ, Isj, 1^ as well as the internal packing 
arrangements of semiflexible polymer torii [3, ,^ ,^] much 
remains that is not understood. Our intention in this 
paper is to highlight the possibilities that arise when one 
extends the problem of semiflexible polymer collapse in 
poor solvent to a copolymer. Inspired by the fact that 
DNA torii occur as a consequence of a competition be- 
tween solvent and bending energies we study copolymers 
with two block types A and B that are immiscible with 
each other and with the solvent. We begin with an ana- 
lytic treatment of the problem which leads us to predict 
several unexpected stable geometries. Using Brownian 
dynamics simulations we demonstrate that such geome- 
tries can self assemble via collapse of a single semiflexible 
copolymer and that their relative energies are in agree- 
ment with the theory. 



THEORY 

As a starting point for our analytic theory we review 
some basic results for a semi-flexible homopolymer. The 
packing arrangement of a semiflexible polymer that op- 
timizes bending within the torus is certainly not a trivial 
matter and has been treated elsewhere [7| . For simplicity 
we will assume a simple garden hose or circular packing. 
We will also assume that the chain enrolls several turns 
in such a way that it is thick enough to give a well de- 
fined toroidal surface, but it is thin enough so that the 
minor radius of the toroid, r, is much less than the ma- 
jor radius, R. For such a thin toroid, every section of the 
chain has approximately the same curvature. Using two 
theorems of Pappus, the surface area and volume of the 



toroid are: A = iir^Rr and V — I-k^Rt^ . Therefore, the 
free energy of the chain can be written as: 
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(1) 



Here the first term represents the bending energy while 
the second term is the surface energy, L = Nl is the total 
length of the polymer along the chain, N is the number of 
repeating units, I is the monomer length, k is the bend- 
ing constant, and 7 is the surface tension (which is in 
this model the parameter that characterizes the polymer- 
solvent interaction). 

Using the volume constraint V — La^ ( a^ is the cross 
section of the chain ) the minor radius of the toroid can 

= — |-W-^ and equation ^ transforms 



be written as, 
to: 
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If the free energy is minimized with respect to the major 
radius of the toroid, the following dependence is obtained: 
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Here we have used the dimensionless parameters a 
and b=-. 

a 

Consider now a more complex problem: a semi-flexible 
copolymer formed by a sequence of two different immis- 
cible blocks (which we call "block A" and "block B"). 
For simplicity, an equal size for both blocks will be as- 
sumed. If the repeating unit along the chain (that is, the 
sequence of a block A followed by a block B) is called a 
diblock, then, for a copolymer, I will denote the length 
of a diblock while N will denote the total number of di- 
blocks along the chain. The repulsive interaction between 
blocks A and B is assumed to be dominant over all other 
interactions such that it may be modeled by imposing 
the restriction that a block A cannot have lateral con- 
tact with a block B. According to this restriction, if the 
copolymer still forms a toroid in poor solvent conditions, 
the ratio of the toroid perimeter, 2ttR, and the diblock 
length, I, should be necessarily a whole number; that is, 

i? = m^, m = l,2,3, ... (4) 



Here the number m denotes the number of diblocks per 
turn in the toroid. The free energy of the toroid still has 
the same form as in equation |21 and since block types 
A and B will allways have equal surface areas within 
our model we can make the interpretation 7 = LlA±IMl^ 
where 7 is the average surface tension of block types A 
and B. For simplicity we also assume that ka = ^b so 
that only a single k is necessary. We should note how- 
ever, that the value of the radius R is now restricted to 
be chosen from a set of discrete values, given by equation 
0] Hence the free energy itself varies in a discrete way: 
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FIG. 1: Schematic diagram of a polymer in the figure 8 con- 
figuration with n = 2 



The transition from a toroid with m + \ diblocks per 
turn, to a toroid with m diblocks per turn occurs when 
the difference in free energy, (AF)^_,_^ — (AF)^, changes 
from negative to positive. If we define the parameters T 
and /m as in equationsElandElthen this will occur when 
T</™. 
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In general, for tti > 1 , the toroid with m diblocks per 
turn will be stable if the condition |H1 is satisfied. 



/jn-l <^ < frr 



(8) 



When m ^ 1 the discrete jumps in R should be small 
compared to i?, so that we should recover the behavior of 
a homopolymer, given by equation|31 Note that when m 

is big, we can approximate fm ~ ^^ and the transition 

5 
points occur when ^^11 ~ X. Therefore the radius of 

the toroid in the transition, given by the restriction 01 

2 

becomes — = m^T- « (-^1 fosiVs coinciding with 

a 2tt \\/27r/ ' *= 

equation O as expected. 

However, when the parameter T is small (that is, when 
the surface energy dominates over the bending energy) 
the radius of the toroid will reach its lower limit, R = 7^, 
corresponding to m — \ diblocks per turn. This lower 
limit to R also implies a lower limit to the surface area 
within the toroidal geometry. Thus we expect that as 
T becomes still smaller and surface energy even more 
important the polymer must explore alternative confor- 
mations to the simple torus. 

One possible conformation, for small values of the pa- 
rameter T, would be a figure 8 conformation, such as that 
shown in figure „ In this conformation a block A enrolls 
a whole number of turns, n, around a circle of radius p 
in the left side, the subsequent block B enrolls n turns 
around a circle of the same radius in the right side, and 



so on. The free energy of the copolymer in this figure 8 
conformation is approximately equal to the free energy 
of two separated toroids of radius p, but each of them 
with a total chain length -^ instead of L. Therefore, the 
expression for the free energy, AFg, corresponding to the 
figure 8 conformation, would be: 



AFs = -—: + ATT'-iaJLp 



(9) 



Now the restriction of no lateral contact between blocks 
A and B is imposed. Since each block (of length |) enrolls 
n turns (with n a whole number) around the circle of 
radius p, this leads to: 
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(10) 



Lets forget momentarily this last restriction and study 
the associated problem of a homopolymer that somehow 
(for example, by imposing some unspecified restrictions) 
has acquired a figure 8 conformation with free energy 
given by equation |51 The equilibrium radius, p, of both 



sides of the figure 8 would be: p 
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tuting this result in equation IHl we obtain a free energy 
AFg = I (27ra7)^ k'sL's . Comparing with the analogous 

free energy of the toroid conformation of a homopoly- 

4 

mer, AF = | (\/27ra7) ^ k'sL's , we can readily see that 
AFg is always bigger than AF by a factor of 2= . There- 
fore, in the case of a homopolymer, where the value of 
the radius is not restricted, the toroid conformation is 
always more stable than the figure 8 conformation. This 
is not surprising, since there is obviously an increase in 
bending when going from a torus to a figure 8. As noted 
previously, however, in the case of a semi-fiexible copoly- 
mer, for small values of the parameter T, i.e. when the 
surface energy dominates over the bending energy, the 
toroid major radius is not allowed to be smaller than 
i? = 2^, while the radius, p, of both sides of the figure 
8 conformation can still be as small as needed, as long 
as they satisfy the restriction^] Thus for small T it is 
possible for the figure 8 to be more stable than the torus. 



The free energy of a figure 8 conformation with n turns 
per diblock is obtained by substituting equation 1101 into 
equation 13 leading to: 
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A figure 8 conformation would be more stable than a 
toroid conformation if the free energy of the toroid with 
one diblock per turn, (AF)^, is bigger than the free en- 
ergy of a figure 8 conformation, {AFg)^ for some value 
of the number of turns per diblock, n. Using equations^ 
and II II this condition leads to: 
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Note that 51 = and hence there cannot be a transition 
from a toroid conformation to a figure 8 conformation 
with just one turn per diblock, n — I. However, for 
values of the parameter T smaller than 32 ~ 0.0195, a 
transition to a figure 8 conformation with two turns per 
diblock would be favourable. 

Finally, transitions between figure 8 conformations 
with different values of n can also be studied by analyz- 
ing the difference between their free energies, (Ai^8)„ ~ 
{AFs)j^_^-^. When this difference becomes bigger than 
zero, the conformation with n + 1 turns per diblock be- 
comes more stable than the conformation with just n 
turns per diblock. It is easy to see that 



{AFs)„ - {AFs)„^, > 
1 \ 1 
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Summarizing the previous discussion the following 
regimes are expected: 

(i) In the case in which inequality |H1 is valid, a toroidal 
conformation with m diblocks per turn (for m > 1). The 
radius of the toroid would be given by equation 01 

(ii) If (72 ^ T ^ /i J then a toroidal conformation with 
one diblock per turn would be expected. The toroid ra- 
dius would be given by i? = ■^. 

(iii) li h2 < T < g2, then a figure 8 conformation with 
one turn per block is expected. 

(iv) Finally, a figure 8 conformation with n turns per 
diblock may be expected if hn+i < T < hn- 



BROWNIAN DYNAMICS SIMULATIONS 

As a compliment to the analytic treatment above we 
have conducted Brownian dynamics simulations with the 



intention of confirming that structures such as the figure 
8 and toroids with to > 1 are realistic polymer mor- 
phologies, and that their relative energies conform with 
our theoretical predictions. 

In our simulation model, poor solvent conditions were 
imposed by using a Lennard Jones pair potential between 
like monomers ie Uaa = Ubb = 4e((cr/r)^2 _ (^g.^j.^6^ ^ 
where r is the inter-monomer distance, e is the strength 
of the interaction and is our unit of energy, a is roughly 
equal to the equilibrium inter-bead distance and is there- 
fore related to the parameter a in our analytic work. 
We chose a — 2 for all simulations. For simplicity in 
our analytic work we specified that blocks of A type 
monomers and B type monomers cannot have lateral 
contact. Of course, such a rigid constraint cannot oc- 
cur for a real polymer and a more realistic situation is 
that monomers A and B experience a repulsive interac- 
tion which offsets their solvent induced attraction. We 
model this in our simulation by specifying the pairwise 
excluded volume interaction between monomers A and B 
as, C/As=4e(a/2r)i2. 

In order to account for the stiffness of the chain we used 
the potential Ubend — ^k'cosO where 9 is the angle be- 
tween successive bonds and k' is the bending constant, k' 
is related to the quantity k used in analytic work through 
the relation k' = ku. 

Bonds between successive monomers were modeled 
using the potential Ubond — kTd~^i*{r/{dj) which is 
equivalent to a Gaussian chain entropic spring potential 
3kTr'^ / {2cP) at short extensions but increases rapidly for 
extensions close to or beyond the bond length d. For 
simplicity, we chose a bond-length equal to the bead di- 
ameter such that a Ki a — d — 2. We calculated par- 
ticle motion by integrating the Langevin equation with 
a time-step of At = 0.002 a temperature of ksT = 1, a 
frictional coefficient C = 1 ^'^d random fluctuations f{t) 
given by a Gaussian distribution with mean and vari- 
ance \j2kBT At. All simulations were carried out for a 
total of 25 X 10^ time-steps. 

Since our intention was to observe striped toroidal 
structures with different numbers of blocks per turn, ra 
we initially chose to simulate a polymer with a large num- 
ber of blocks N = 20. The total chain length was 400 
monomers, which gives L = 800 and I = LjN = 40. 
Stiffness was chosen to be relatively large k! — 15 so 
as to avoid globular configurations which are not ac- 
counted for in our analytic model. We conducted 40 
replicate runs in which a polymer with the above param- 
eters was quenched from an equilibrated configuration 
in good solvent conditions to poor solvent via the im- 
position of inter-monomer LJ potentials. We found that 
after an initial period of collapse, a steady state configu- 
ration was reached which showed no change over at least 
the final half of the simulation time. A wide variety of 
such configurations were observed, including those with 
TO = 2,3,4 as well as some previously unexpected mor- 
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FIG. 2: Steady-state collapsed configurations of semiflexible 
copolymer where A'^ = 20, L = 400, k, = 15. Snapshots shown 
illustrate the diversity of collapsed states and are labeled al- 
phabetically. The relative internal energies of these states, 
t/(e) and their frequency of occurrence / as a percentage of 
the 40 replicate simulation runs are as follows. A: U — 0±10, 
f = 20%, B: f/ = 115 ± 20, / = 10%, C: f/ = 30 ± 20, 
/ = 20%, D: f/ = 250 ± 10, / = 12.5%, E: U ^ 360,/ = 2.5%, 
F: U = 300 ± 40, / = 27.5% . Overall values for U and 
associated errors are calculated from the set of all replicates 
that collapsed to a particular configuration. Values for each 
of these simulations in turn were taken over the final 1 x 10^ 
time steps of the simulation. 



phologies. Since entropic effects are negligible in our poor 
solvent system the free energies of each of these config- 
urations are roughly equal to their internal energies U. 
Under this assumption we found that the m = 2 config- 
uration was the most stable of those observed, with all 
other states including toroids with m = 3 and 4 all being 



higher in energy. 

Close examination of the configurations in figure „ re- 
veals several points of particular note. Firstly, we see 
that the repulsion between A and B type monomers ef- 
fectively results in a complete segregation of these blocks. 
The end result is almost identical to the constraint used 
in our analytic calculations that A and B blocks cannot 
have lateral contact. What is not predicted in our an- 
alytic theory is that this segregation of A and B blocks 
leads to the formation of a diverse range of metastable 
structures. This occurs because A and B type monomers 
cannot slide against one another, which is an important 
equilibration mechanism for homopolymers. As a conse- 
quence, we find that structures with m = 3 and m — A 
cannot readily transform to the more stable m — 2 type 
torus. Instead they remain as metastable states. In ad- 
dition to the striped toroidal states we also observe other 
interesting metastable structures such as a torus within 
a torus (figure bC) or a torus at the end of a rod (figure „ 
E) or a striped rod (figurCnF). In all cases it seems likely 
that the A B segregation prevents a transition to a state 
of lower free energy. 

Although the simulation parameter e is clearly related 
to the surface tension 7, this relationship is very difficult 
to obtain from a first principles calculation. Instead, we 
can calibrate 7 for our simulation system by using the 
internal energies of the striped toroidal states. The rela- 
tive energies of these states should be given by equation 
Elso by fitting this equation to values of U and allowing 
for an unknown additive constant, we can obtain a rough 
estimate for 7. The value of 7 so obtained was 7 « 0.4. 
By substituting this value of 7 back into equation \^ we 
find that the energy of the m = 1 torus for our simu- 
lated polymer is predicted to be C/ « 800e. Such a high 
predicted free energy is consistent with the fact that we 
never observed the formation of to = 1 toroids for the 
TV = 20 polymer. Our value of 7 can also be used to 
predict an appropriate set of parameters for the observa- 
tion of figure 8 structures. In making such a prediction 
we keep the total polymer length L constant, such that 
I = L/N, and vary TV. Keeping all other parameters 
constant we predict that figure 8 configurations should 
be observed for TV < 6. Taking a conservative approach 
to this prediction we simulated a polymer with TV = 2 
using the same methods as were used for the TV := 20 
case. 

In figure „ we present snapshots of steady state col- 
lapsed configurations of the N ~ 2 polymer. The most 
stable observed configuration in this case is the figure 8, 
however as with the TV = 20 polymer many metastable 
configurations were also observed. Indeed the most com- 
monly observed shape was that of a single torus attached 
to a rod. Other metastable states included complete rods 
as well as an ?n = 1 toroid and a strange configuration 
in which a toroid exists between two ends of a bent rod. 
Although this last configuration is rare it serves to illus- 
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FIG. 3: Steady-state collapsed configurations of semiflexible 
copolymer where N — 2, L — 800, k = 15 . Snapshots shown 
illustrate the diversity of collapsed states and are labeled al- 
phabetically. The relative internal energies of these states, 
(7(e) and their frequency of occurrence / as a percentage of 
the 40 replicate simulation runs are as follows. A: t/ = 70±10, 
/ = 37.5%, B: f/ = 0±10, / = 15%, C: U = 440±10, / = 5%, 
D: (7 = 140, / = 2.5%. Rods [7 = 220 ± 10,/ = 35% are not 
shown. See figure „ for details on calculation of averages and 
errors 



trate the diversity of metastable states which occur when 
A and B type blocks repel one another. Provided the A 
B repulsion is sufficiently strong we would expect that 
many of the the observed metastable states should be 
quite long lived. 



CONCLUDING REMARKS 



tures in poor solvent. Even the very simple alternating A 
B block case discussed here forms a range of torii with dis- 
crete sizes as well as the unusual figure 8 structure. Aside 
from this variety of structures one of the key features of 
these copolymer collapsed states is their well defined size 
which must be an integer multiple of the block length. 
Although such structures are yet to be realized we ex- 
pect that they could have great utility as tools in single 
molecule experiments, where for example a molecule is 
pulled through a precisely defined pore ( eg torus ) llQj. 
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Semiflexible copolymers exhibit the potential for for- 
mation of a wide variety of single molecule compact struc- 



